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FREQUENCY  STABILITY  SPECIFICATION  AND  MEASUREMENT: 
HIGH  FREQUENCY  AND  MICROWAVE  SIGNALS* 

John  H.    Shoaf,    Donald  Halford,    and  A.    S.    Risley 

Atomic  Frequency  and  Time  Standards  Program  Areas 
National  Bureau  of  Standards 
Boulder,    Colorado     80302     USA 

This  report  gives  concise  definitions  for  specifying 
frequency  stability  for  measurements  in  the  frequency 
domain  and  time  domain.      Standards  of  terminology  and 
of  measurement  techniques  are  recommended.     Measure- 
ment systems  in  the  high  frequency  and  microwave  regions 
are  described  in  adequate  detail  so  that  the  systems  may 
be  duplicated. 

Key  Words:     Allan  variance;  Frequency  stability  measure- 
ments; Measurement  system  description; 
Phase  noise;    Spectral  density;  Stability 
definitions;  Terminology  standards. 


I.  INTRODUCTION  AND  BACKGROUND 

At  the  beginning  of  FY- 71  the  Department  of  Defense  Joint 
Services  Calibration  Coordinating  Group  (DoD/CCG),    J.    L.    Hayes 
(Chairman),    Metrology  Engineering  Center ,    Pomona,    California, 
contracted  with  the  National  Bureau  of  Standards   (NBS)  to 
write  a  paper  pertaining  to  the  specification  and  measurement 
of  frequency  stability.      The  project  was  under  the  jurisdiction  of  the 
DoD/CCG  Time  and  Frequency  Working  Group.      The  current  members 
of  this  group  are  Peter  Strucker   (Chairman),    Metrology  Engineering 
Center,    Pomona,    California;  J.    M.    Rivamonte,    U.    S.    Army  Metrology 
and  Calibration  Center,    Redstone  Arsenal,    Alabama;  and  E.    L. 
Kirkpatrick,    Aerospace  Guidance  and  Metrology  Center,    Newark  Air 
Force  Station,    Newark,    Ohio. 


*  Contribution  of  the  National  Bureau  of  Standards,    not  subject  to 
copyright. 


The  first  year's  work  (FY-71)  culminated  as  NBS  Report  9794  [1]. 
The  second  year's  work  (FY -72)  extended  the  frequency  range  to  include 
a  working  measurement  system  at  X-band.  The  results  of  work  done  in 
both  FY-71  and  FY -72  are  documented  in  this  paper. 

The  purpose  of  the  project  is  to  establish  recommended  standards 
of  terminology  and  measurement  techniques  for  frequency  stability. 
Emphasis  is  placed  on  details  of  useful  working  systems   (apparatus) 
that  could  be  duplicated  by  others  in  the  field  of  frequency  stability 
measurements.     Uniformity  of  data  presentation  is  stressed  in  order  to 
facilitate  interpretation  of  stability  specifications  and  to  enable  one  to 
communicate  and  compare  experimental  results  more  readily.      An 
authoritative  paper    [2]  by  the  Institute  of  Electrical  and  Electronics 
Engineers   (IEEE)  Subcommittee  on  Frequency  Stability  is  used  exten- 
sively as  the  prime  reference,  in  the  preparation  of  this  report.      This 
report  presents  the  description  and  performance  of  frequency  stability 
measurement  systems  capable  of  precise  measurements  on  state-of-the- 
art  sources.      Sufficient  theory  of  the  measurement  systems  is  given  to 
enable  a  person  to  readily  understand  the  principle  of  operation. 

2.     TERMINOLOGY  FOR  SPECIFICATION 
OF  FREQUENCY  STABILITY 

The  term  frequency  stability  encompasses  the  concepts  of  random 
noise,    intended  and  incidental  modulation,    and  any  other  fluctuations  of 
the  output  frequency  of  a  device.     In  this   report  we  are  mainly  (but  not 
totally)  concerned  with  random  fluctuations  corresponding  to  Fourier 
frequencies  in  the    10     to   10     hertz  range.      The  measurement  of  fre- 
quency stability  can  be  accomplished  in  both  the  frequency  domain  (e.  g.  , 
spectrum  analysis)  and  the  time  domain  (e.  g.  ,    gated  frequency  counter). 


In  the  aforementioned  manuscripts    [1],  [2]  the  authors  chose  to  use 

two  independent  definitions,    each  related  to  different  useful  methods  of 

measurement.      (See  Appendix  A  for  a  Glossary  of  Symbols.  ) 

The  frequency  domain  definition  of  frequency  stability  is  the 

one-sided  spectral  density  of  the  fractional  frequency  fluctuations,    S    (f), 

where   y   =    glV  V  •     The  fractional  frequency  fluctuation  spectral  density 

S    (f)    is  not  to  be  confused  with  the  radio  frequency  power  spectral 

density   S  / (v),  nor  with   S       (v) ,    which  are  not  good  primary  meas- 

\'RFP  oV  

jares  of  frequency  stability  [2],      (There    is  some  discussion  of  this  in 

Appendix  D.  )     Phase  noise  spectral  density  plots  [i.  e.  ,    S     ,  (f)  versus  f] 

00 

are  a  common  alternative  method  of  data  presentation.     The  spectral 
density  of  phase  fluctuations  is  related  to   S   (f)   by 

s*Jf>   =   (f2/f2)s  (f)    •  (1) 

60  o  y 


The  time  domain  definition  of  frequency  stability  uses  the  type  of 
sample  variance  called  the  Allan  variance    [3]  of  y: 

<*>T,r,fh»    .  /^f  ^.JL^^    .  m 

\  n=l  k=l 


The  particular  Allan  variance  with  N  =  2   and   T  =  T  is  found  to  be 
especially  useful  in  practice.     It  is  denoted  by: 

/(yk+1  -  yk> 

</(T)     =      <0y(N  =  2,T=T,r,fh)>    :      (  , 

The  bar  over  the  y  indicates  that  y  has  been  averaged  over  a  specified 
time  interval  T.  The  angular  brackets  indicate  an  average  of  the  quan- 
tity over  time.      (See  example  in  Appendix  F.  ) 


In  the  time  domain  we  are  concerned  with  the  measurement 

2         ^ 
of  Allan  variances  at  different  time  intervals.  Plots  of    [ex    (T)]2 

versus    T    ("sigma  versus  tau")  on  a  log-log  scale  are  commonly  used 
for  data  presentation.     A  convenient  chart  which  enables  one  to  trans- 
late from  frequency  domain  measures  to  time  domain  measures   (and 
often  conversely)  is  found  in  Appendix  B.     An  example  of  this  transla- 
tion is  given  in  Appendix  C. 

Script   ?£(f)   is  a  frequency  domain  measure  of  phase  fluctuations 
(noise,   instability,   modulation)  used  at  NBS.     Script   s£(f)   is  defined  as 
the  ratio  of  the  power  in  one  phase  noise  sideband,    referred  to  the  input 
carrier  frequency,    on  a  per  hertz  of  bandwidth  spectral  density  basis, 
to  the  total  signal  power,    at  Fourier  frequency  f  from  the  carrier,   per 
one  device.     (See  Appendix  D.  ) 

Power  density  (one  phase  modulation  sideband) 
(  '     =  Power  (total  signal)  *  (   ) 


For  small    60, 


S       (f)    =    2£(f)  .  (5) 


A  practical  system  for  the  measurement  of  Script  £(i)    or  sfiJf)   wi11  be 
described  in  detail  later. 

It  seems  appropriate  here  to  discuss  briefly  the  types  of  noise 
that  affect  the  output  frequency  of  all  known  signal  sources.     The  noises 
can  be  characterized  by  their  frequency  dependence. 

Common  Name 

White  Phase  Noise 

Flicker  Phase  Noise 

White  Frequency  Noise 
(Random  walk  of  phase) 

Flicker  Frequency  Noise 


S60(f) 

Vf)  =   f  So0(f) 

f° 

f+2 

f"1 

f+1 

f"2 

f° 

f"3 

f"1 

3.     OPERATIONAL  SYSTEMS  FOR  MEASUREMENT 

OF  FREQUENCY  STABILITY  AT  NBS 

(HIGH  FREQUENCY  REGION) 

John  H.    Shoaf  and  Donald  Halford 

Most  conventional  systems  for  measurement  of  frequency  stability 
until  recently  have  utilized  time  domain  techniques  primarily.     As 
indicated  previously,   in  order  to  have  a  comprehensive  and  sufficient 
measure  of  frequency  stability,   it  is  preferred  that  the  measurements 
involve  both  frequency  and  time  domain  techniques. 

Fortunately,   frequency  domain  and  time  domain  methods  for 
measuring  frequency  stability  require  similar  apparatus  except  that: 
to  make  measurements  in  the  frequency  domain  you  must  have  a 
frequency  window  (spectrum  analyzer)  following  the  detector;  for  the 
time  domain  you  must  have  a  time  window  (gated  counter)  following  the 
detector. 

It  was  the  introduction  of  good  double -balanced  mixers  that  per- 
mitted measurement  of  frequency  stability  by  improved  techniques  [4] -[7], 
[11].     The  double -balanced  mixer,    considered  as  a  phase  sensitive  de- 
tector,  makes  possible  meaningful  frequency  stability  measurements  of 
high-quality  signal  sources  in  both  the  frequency  domain  and  the  time 
domain.     The  results  are  quantitative  and  may  be  obtained  from  a 
measurement  system  which  is  reasonable  in  cost. 

The  frequency  stability  measurement  systems  described  below 
have  been  used  at  NBS  since   1967.     The  functional  block  diagrams  in 

figures   1,    7,  9,  and  10  are  referred  to  in  the  detailed  descriptions  of  the 

3  9 

particular  systems.       The  carrier  frequency  range  10     Hz  to  10     Hz  is 

easily  covered  with  these  techniques. 


3.  1.     Frequency  Domain  Measurements 
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FIGURE  1;  TYPICAL  FREQUENCY  DOMAIN  MEASUREMENT  OF  FREQUENCY  STABILITY 
(PHASE  SENSITIVE  MODE) 


OSCILLATOR  UNDER  TEST 

REFERENCE  OSCILLATOR  OF  HIGH  SPECTRAL  PURITY  (VARACTOR  TUNING) 

ADJUSTABLE  ATTENUATOR  (TYPICALLY    lOdB) 

DOUBLE-BALANCED  MIXER  (SCHOTTKY  BARRIER  DIODES) 

NBS  LOW-NOISE  dc  AMPLIFIER  (SEE  FIGURE  3) 

OPERATIONAL  AMPLIFIER  (SEE  FIGURE  2) 

BATTERY  BIAS  BOX  (SEE  FIGURE  5) 

SPECTRUM  ANALYZER  (LOW  FREQUENCY,  HIGH  RESOLUTION) 


(Figures   1,    7,    and  10  show  circuit  values  appropriate  for 
measuring  high  quality  5 -MHz  quartz  crystal  oscillators.  ) 


Figure  1  illustrates  the  measurement  system  typically  used  at 
NBS  for  frequency  domain  measurements.     It  should  be  noted  that  time 
domain  data  can  also  be  obtained  simultaneously,    although  usually  this 
system  is  used  only  for  frequency  domain  measurements.      (For  time 
domain  measurements  it  is  often  more  convenient  to  use  a  slightly 


modified  measurement  setup  to  be  described  later.  )     In  this  frequency- 
domain  setup  the  oscillator  under  test  is  fed  into  one  side  of  a  low-noise 
double-balanced  mixer  which  utilizes  Schottky  barrier  diodes.     The 
reference  oscillator  is  fed  into  the  other  side  of  the  mixer  through  an 
attenuator,    typically  10  dB.      The  mixer  acts  as  a  phase  sensitive 
detector  so  that  when  the  two  signals  are  identical  in  frequency  and  are 
in  phase  quadrature  the  output  is  approximately  zero  volts  dc.     When 
this  output  is   sent  back  to  the  reference  oscillator  via  the  varactor 
tuning,    phase  lock  is  achieved.  As  with  any  feedback  system,    care  must 
be  observed  to  avoid  instabilities,    e.  g.  ,    servo  loop  oscillations.      The 
phase-lock  loop  contains  proper  termination  at  the  output  of  the  mixer 
followed  by  operational  amplifiers  with  adjustable  gain.      The  time 
constant  of  the  loop  may  be  adjusted  as  needed  by  varying  the  amplifier 
gain  within  the  loop  and  by  use  of  the  RC  filter  (fig.    4)  indicated  in  the 
diagram.      Finally,    a  battery  bias  box  is  included  at  the  varactor  input 
in  order  to  operate  in  a  suitably  linear  portion  of  the  varactor' s 
frequency  versus  voltage  curve. 

A  very  loose  phase-lock  loop  is  indicated  inasmuch  as  the  voltage 
varies  as  phase  (in  short  term),    and  in  this  frequency  domain  measure- 
ment we  are  observing  the  small  phase    variations  directly.     By  the 
phrase  very  loose  phase-lock  loop,    we  mean  that  the  bandwidth  of  the 
servo  response  is  small  compared  to  the  lowest  frequency  f  at  which 
we  -wish  to  measure  (i.e.,   the  response  time  is  very  slow).     Operational 
amplifiers  are  arranged  in  a  circuit  as  shown  in  figure  2  for  convenience 
of  adjusting  the  gain  and  for  self-contained  battery  supply  voltage. 
Special  NBS  low-noise  dc  amplifiers  used  in  certain  precision  measure- 
ments are  shown  in  figure  3.     At  NBS  we  have  arranged  in  a  small 
chassis  the  adjustable  RC  or  CR  filters  utilizing  low-noise  components, 
with  rotary  switches  for  various  combinations  of   R   and   C    (see  fig.    4). 
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FIGURE  2:  STEPPED  GAIN  OPERATIONAL  AMPLIFIER 
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The  battery  bias  box  is  arranged  with  vernier  as  shown  in  figure  5, 
facilitating  fine  frequency  adjustments  via  the  varactor  frequency  adjust- 
ment in  one  oscillator.     A  wave  analyzer  is  used  to  obtain  the  noise  plot 
information  relevant  to  stability  (frequency  domain).      The  phase  noise 
sideband  levels  are  read  out  in  rms  volts  on  the  analyzer  set  to  certain 
chosen  values  of  frequency,    f.      For  typical  high  quality  signal  sources 
this  corresponds  to  measuring  only  those  phase  noise  sidebands  which 
are  separated  from  the  carrier  by  the  various   f   intervals  chosen. 
Script   =£(f)   may  be  calculated  with  the  assumption  that  both  sources 
contribute  equally;    however,    if  one  source  -were  the  major  contributor, 
then  the  noise  of  that  source  would  be  no  worse  than  3  dB  greater  than 
the  value  of  Script    s£(f)   so  calculated.     A  typical  plot  is  shown  in  figure  6. 
A  sample  calculation  may  be  found  in  Appendix  E. 
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3.  2.      Time  Domain  Measurements 


In  figure  7  a  measurement  system  typically  used  at  NBS  for 
stability  measurements  in  the  time  domain  is   shown.     It  will  be  noted 
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FIGURE. 7:  TYPICAL  TIME  DOMAIN  MEASUREMENT  OF  FREQUENCY  STABILITY 
(FREQUENCY  SENSITIVE  MODE) 

ITEMS  (1)  THROUGH  (8)  SAME  AS  FIGURE  1 

(9)  OPERATIONAL  AMPLIFIER 

(10)  STRIP  CHART  RECORDER  FOR  QUALITATIVE  OBSERVATION 

(11)  VOLTAGE-TO-FREQUENCY  CONVERTER 

(12)  FREQUENCY  COUNTER  WITH  LOW  DEAD  TIME 

(13)  DIGITAL  RECORDER  WITH  FAST  RECORDING  SPEED 
(INHIBIT  TIME  COMPATIBLE  WITH  COUNTER  DEAD  TIME) 

(14)  COMPUTER  (OPTIONAL  METHOD  OF  DATA  ANALYSIS) 

that  the  principle  of  operation  is   similar  to  that  used  in  the  frequency 
domain  measurement  wherein  the  reference  oscillator  is  locked  to  the 
test  oscillator.     However,    for  the  time  domain  measurement  we  use  a 
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very  tight  phase-lock  loop  and  the  correction  voltage  at  the  oscillator 
varies  as  frequency.     By  the  phrase  very  tight  phase-lock  loop,    we  mean 
that  the  bandwidth  of  the  servo  response  is  relatively  large  (i.  e.  ,    the 
response  time  is  much  smaller  than  the  smallest  time  interval    T  at 
which  we  -wish  to  measure).       Caution- -there  are  problems  'which  are 
potentially  present  in  tight  phase-lock  loop  systems.      This  is  a  very 
convenient  setup  for  observing  frequency  fluctuations  in  longer  term. 
However,    with  the  time  constant  appropriately  adjusted  and  the  means 
for  taking  sufficiently  fast  samples  the  system  is  readily  used  for  short 
term  measurements,    as  well  as  for  the  longer  term  measurements,    in 
the  time  domain.      For  qualitative  observations  any  suitable  oscilloscope 
or  strip  chart  recorder  may  be  used.      For  quantitative  measurements 
the  system  at  NBS  utilizes  a  voltage-to-frequency  converter,    a  frequency 
counter,    and  a  printer  capable  of  recording  rapid  samples  of  data  with 
very  short  dead  time.      The  data  are  analyzed  typically  by  computer  via 
a  program  designed  to  compute  the  appropriate  Allan  variance    [2],  [3]. 
In  our  computer  program  log    cr   versus  log    T  along  with  the  associated 
confidence  in  the    a   are  automatically  plotted  on  microfilm.      For  small 
batches  of  data  a  desk  calculator  could  be  used  and  the  computer  analysis 
would  not  be  necessary.      An  example  of  a  specific  Allan  variance  compu- 
tation is  shown  in  Appendix  F.      A  typical  plot  of  log   a   versus  log    T  is 
shown  in  figure  8.      The  dashed  lines  indicate  the  slopes  which  are 
characteristic  of  the  types  of  noise  indicated. 
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The  convenience  of  obtaining  time  domain  data  has  been  greatly- 
enhanced  by  utilizing  recently  developed  counters  [8]  which  are  program- 
able  to  automatically  compute    a    versus    T.     A  block  diagram  (fig.    9) 
shows  the  measurement  setup  using  a  computing  counter  programed  for 
(T    (t)  .      The  program  -which  was  used  is  given  in  Appendix  G. 
Certain  limitations  of  deadtime  are  inherent  in  the  use  of  this  time 
domain  method.     However,    in  general  (except  for  very  short    T), 
frequency  stability  in  the  time  domain  may  be  measured  quickly  and 
accurately  using  a  computing  counter. 
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MIXER 

FILTER 

COMPUTING 
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PRINTER 

i 

SIGNAL 
SOURCE 

FIGURE  9:  FREQUENCY  STABILITY  MEASUREMENT  UTILIZING  A  COMPUTER  COUNTER 
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3.  3.     Differential  Phase  Noise  Measurements 

An  additional  useful  system  illustrated  in  figure    10  is  used  for 
differential  phase  noise  measurements  of  various  discrete  components 
which  are  frequently  used  in  stability  measurement  systems.     In  this 
system  only  one  frequency  source  is  used.     Its  output  is   split  so  that 


part  of  the  signal  passes  through  the  component  to  be  tested. 
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FIGURE  10:  DIFFERENTIAL  PHASE  NOISE  MEASUREMENT 
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ITEMS  (1)  THROUGH  (14)  SAME  AS  FIGURES  1  AND  7 

(15)  ANY  DEVICE  OR  COMPONENT  UPON  WHICH  NOISE  MEASUREMENTS  ARE  DESIRED 
(AMPLIFIERS,  FILTERS,  CAPACITORS,  CABLES,  PADS.  ETC.) 

(16)  NBS  ADJUSTABLE  PHASE  SHIFTER,  5MHz  (SEE  FIGURE  11) 

The  signal  is  adjusted  via  a  phase  shifter   (fig.    11)   so  that  it  is  in  phase 
quadrature  with  the  other  part  of  the  original  signal  and  is  down- 
converted  in  the  Schottky  barrier  diode  mixer  as  described  in  the  other 
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FIGURE  11:  ADJUSTABLE  PHASE  SHIFTER  (5MHz  DELAY  LINE) 

RG174/U  CABLE  WAS  USED  FOR  EACH  SEGMENT  OF  PHASE  SHIFT 
CALCULATED  AT  ~ 1 0  cm  PER  DEGREE  AT  5  MHz. 
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systems.     The  switchable  50f2  load  in  figure  II  is  not  essential  but  is 
included  for  convenience.     A  low-pass  filter  is  included  before  the 
signal  is  amplified  in  special  low-noise,   low-level  dc  amplifiers  and 
observed  on  the  spectrum  analyzer.     Script   :£(f)   values  are  calculated 
at  various  frequency  values,    f,    and  plotted.     A  sample  calculation  is 
shown  in  Appendix   E. 

The  measurement  system  noise  level  (e.  g.  ,    see  fig.    6)  is  easily 
evaluated.     Using  the  differential  phase  noise  measurement  system 
shown  in  figure  10,  let  the  "device  under  test"  be  a  short  length  of  coaxial 
cable     (which  is  itself  not  a  source  of  noise).      The  small  amount  of  noise 
observed  on  the  spectrum  analyzer  represents  the  system  noise,    mainly 
due  to  the  mixer  (4)  or  the  first  amplifier   (5).      The  calculation  of  the 
system  noise  is  then  the  procedure  given  in  Appendix  E. 
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4.     OPERATIONAL  SYSTEMS  FOR  MEASUREMENT 
OF  FREQUENCY  STABILITY  AT  NBS 
(MICROWAVE  REGION) 

John  H.    Shoaf  and  A.    S.    Risley 

Thorough  investigation  of  stability  measurement  techniques  in 
the  X-band  region  revealed  that  a  method  different  from  that  described 
above  for  HF  measurements  -was  desirable  in  most  cases.     The  recom- 
mended measurement  system  is  described  here  in  detail.      Other 
techniques  of  stability  measurements  in  X-band  will  be  discussed  in 
less  detail  later.     The  recommended  system*  is  a  single -oscillator 
system  as  shown  in  the  photograph  (fig.    12)     and  in  the  block  diagram 
(fig.    13).     Prime  references  are  papers  by  Ashley  et  al.  [9]  and 
Ondria   [10]. 

4.  1.  Discussion  of  the  Measurement  System 

The  single-oscillator  frequency  stability  measurement  system  is 
basically  a  frequency  modulation  (FM)  demodulator.      That  is,    it  can 
retrieve  from  the  modulated  carrier  the  signal  with  which  the  carrier 
was  originally  frequency  modulated. 

An  important  consideration  when  making  measurements  is  that  of 
maintaining  the  quadrature  condition- -a  90°  average  phase  difference 
between  the  signals  in  the  reference  channel  and  the  signal  channel  as 
seen  at  the  mixer.     Unfortunately,    there  is  a  fairly  high  probability  that 
during  the  course  of  a  measurement  the  average  phase  difference  will 
fluctuate  a  few  degrees  about  the  desired  90°   setting.     Therefore,    it  is 
recommended  that  an  occasional  check  of  the  quadrature  condition  be 


*  The  recommended  system  is  discussed  here  as  a  frequency  domain 
measurement.     However,    time  domain  measurements  can  also  be  made. 


made.      (In  a  two- oscillator  system  of  measurement  discussed  later  the 
quadrature  condition- -in  long  term- -is  established  and  maintained  by 
phase-locking  one  source  to  the  other.     A  similar  procedure  could  be 
used  here,    but  we  consider  it  to  be  unnecessary  in  practice.  ) 

In  practice,    there  is  a  low  frequency  limit  to  the  usefulness  of 
this  method  for  the  measurement  of  FM  noise.     We  have  seen  limiting 
values  of   f  'ranging  from  as  high  as  500  Hz  to  as  low  as  2.  5  Hz. 

The  single-oscillator  system  and  the  two-oscillator  system  each 
has  an  upper  frequency  limit;  i.  e.  ,  a  value  of  f  above  which  frequency 
stability  measurements  cannot  meaningfully  be  made.  For  the  single- 
oscillator  system  this  upper  limit  is    f  «    W    ,    where    W      is  the  3-dB 

c  c 

resonance  line-width  of  the  loaded  discriminator  cavity.     In  the  NBS 
single-oscillator  system,    measurements  were  made  at  values  of   f   as 
high  as    100  kHz. 
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4.  2.      Description  of  the  Measurement  System 
A  description  of  the  single -oscillator  measuring  system  may 
readily  be  followed  by  referring  to  the  block  diagram  (fig.    13).      The 
X-band  source  under  test  is  connected  at  the  left-hand  side  of  the  sys- 
tem.    The  signal  passes  through  an  isolator    {l}  and  variable  attenuator 
{2}  before  it  is  split  via  a  3-dB  directional  coupler    {3}.      (It  should  be 
noted  that  isolators    {l},  {4},  {7},  {ll},    and   114}  are  used  at  several 
points  throughout  the  system  as  a  means  of  preventing  any  serious 
reflections  which  might  otherwise  exist.  )     Part  of  the  signal  enters  the 
reference  channel  (upper  arm),    passing  through  a  phase  shifter    {6}  via 
a  variable  attenuator    {5}  and  eventually  through  a  90°  twist  into  the 
balanced  mixer    {8}.      The  other  part  of  the  signal  enters  the  signal 
channel  (lower  arm),    passing  through  a  three-port  circulator    {9}  con- 
nected to  a  discriminator  cavity   {10}  at  one  port.     A  variable  attenuator 
[12}  is  also  in  the  signal  channel  before  the  signal  reaches  the  balanced 
mixer.     The  output  of  the  mixer  goes  to  either  of  two  spectrum  analyzers 
{17}.      A  10-dB  directional  coupler    {l3}is  utilized  in  the  signal  channel 
to  facilitate  detection  of  resonance  tuning  of  the  cavity.     This  is  observed 
via  a  detector    {15}   with  a  dc  voltmeter  readout    {16}. 

The  only  component  in  the  system  which  is  not  readily  available 
commercially  as  a  stock  item  is  the  discriminator  cavity    [10].     It  is  a 
TE  right  circular  one-port  cavity.      The  coarse  tuning  is  accomplished 

by  means  of  a  movable  end  wall.      The  fine  tuner  is  a  small  diameter  rod 
■which  can  be  moved  coaxially  in  the  cavity.      The  diameter  of  the  rod 
should  be  such  that  the  cavity  frequency  changes  no  more  than  1.  5  MHz 
for  0.05  inch  (~  1  millimeter)  change  in  depth  of  insertion.      At  any 
desired  frequency  the  coupling  of  the  cavity  should  be  such  that  the 
absorption  is  very  nearly  complete.      For  further  discussion  of  coupling 
see  reference   [4].      The  cavity  Q  needs  to  be  high  enough  for  good 
sensitivity  but  adequately  low  for  sufficient  bandwidth.      The  cavity  used 
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in  the  frequency  stability  measurement  system  described  here  has  an 
unloaded  Q  of  approximately  20,  000.  Additional  details  are  available 
upon  request  from  the  authors. 

4.  3.     Calibration  Procedure 

Initial  calibration  of  the  measurement  system  is  necessary  in 
order  to  assign  an  absolute  scale  to  the  stability  measurements.      To 
facilitate  calibration,    a  sinusoi  dally -modulated  X-band  source  is  used  to 
drive  the  system.      The  frequency-modulated  signal  is  observed  on  an  RF 
power  spectral  density  analyzer  and  the  modulation  level  is  adjusted  to  a 
value  sufficient  to  completely  suppress  the  X-band  carrier.      For  sinus- 
oidal modulation,    the  first  carrier  null  corresponds  to  a  modulation 
index  of  2.  4.     Modulation  at  5  kHz  was  found  to  be  convenient  to  use 
because  of  the  particular  dispersion  and  bandwidth  settings  which  -were 
available  on  the  particular  spectrum  analyzer  used  to  display  the  carrier 
suppression.      The  detailed  procedure  for  obtaining  the  calibration  factor 
follows. 

(a)  With  the  discriminator  cavity   {10}    far  off  resonance,    set  the 
level  of  the  X-band  signal  (as  determined  with  a  dc  voltmeter    {16}  at  the 
detector    {l5})  to  a  convenient  value  and  record  the  value.      The  first 
variable  attenuator    {2}  should  be  used  for  this  adjustment.     Any  con- 
venient level  may  be  chosen  provided  an  equal  amount  of  power  also  will 
be  available  from  the  signal  source  which  is  to  be  evaluated. 

(b)  Adjust  the  cavity  to  resonance.     The  dc  voltmeter  at  the 
detector  is  used  to  determine  resonance.     Place  the  dc  voltmeter  at  the 
output  of  the  mixer    (8}  and  adjust  the  phase  shifter    {6}  in  the  reference 
channel  until  the  dc  output  at  the  mixer  is  zero  (phase  quadrature). 
Remove  the  dc  voltmeter  and  connect  the  output  of  the  mixer  to  a  spectrum 
analyzer    {17}  tuned  to  the  modulation  frequency,    5  kHz.     Record  the  rms 
voltage  reading  (V  )  o-  the  spectrum  analyzer. 
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(c)     It  is  now  possible  to  calculate  the  calibration  factor    K. 

(A  v) 
K    =-7-SS£    ,  (6) 

rms 

where    (A  p)  is  the  rms  frequency  deviation  of  the  carrier  due  to 

rms 

intentional  frequency  modulation.      This  deviation  is  calculated  using  the 

equation 

(Ay)  =    0.707  (Av)        .,  (7) 

rms  peak 

where  {Av)  ,  is  the  product  of  the  modulation  index  with  the  frequency 
of  sinusoidal  modulation,  i.  e.  ,  2.4x5  kHz  =  1Z.  025  kHz.  Therefore  the 
calibration  factor  in  our  case  is 


v        0.  707(12.025  kHz)         8.51  kHz 

K  =  -  =    —  .  (8) 

rms  rms 


4.4    Measurement  Procedure 

The  procedure  for  obtaining  data  for  the  spectral  density  plot  is 
quite  similar  to  the  calibration  procedure  except  that  the  X-band  carrier 
is  not  subjected  to  intentional  modulation. 

(a)  With  the  cavity  far  off  resonance,    set  the  level  at  the  detector 
to  the  same  value  obtained  during  calibration.     Use  the  variable  attenu- 
ator   {2}  for  this  adjustment.      The  other  variable  attenuators    {5},  {12} 
are  set  to  zero. 

(b)  Adjust  the  cavity  to  the  resonant  frequency  of  the  X-band 
source.     Adjust  the  phase  shifter  so  that  the  dc  output  at  the  mixer  is 
zero.     Attach  the  spectrum  analyzer  to  the  output  of  the  mixer  and 
record  rms  voltage  readings   (v  )  for  various  frequency  settings  of 
the  spectrum  analyzer.     A  low  noise  amplifier  may  be  necessary  to 
obtain  useful  readings  at  large  Fourier  frequencies.      A  second  reading 
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(v  )  should  be  taken  at  each  value  of   f   with  the  signal  strongly  atten- 

rms 

uated  in  the  signal  channel.      This  is  to  record  the  residual  additive 
background  noise  not  attributable  to  actual  phase  noise  on  the  carrier. 
This  attenuation  is  accomplished  by  inserting  all  (>  20  dB)  of  the  attenu- 
ation in  the  variable  attenuator    {l2}. 

(c)     In  order  to  calculate  values  of   S      (f)    for  plotting  at  various 
Fourier  frequencies  it  is  convenient  to  make  a  tabulation  of  results.      An 
example  of  some  typical  results  is  given  below  in  Table   1.     The  following 
relations  are  used: 


v 
rms 


V(v'        )2  -  (v"       )2  (9) 

rms  rms 


fiy  =     v  X    K  (10) 

rms  rms 


(6k         )2 

V£)  -  -"it-  <n) 


where   B    is  the  bandwidth  at  which  the  readings  were  made  on  the 
spectrum  analyzer,    and 

S60(f)  =  -J-   ■  (12) 


Values  of   S      (f)    which  were  calculated  this  way  (Table   1)  are  plotted  in 
figure   14. 
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TABLE   1 


Worksheet  for  Calculation  of   S     .(f) 


Calibration  factor  K 


51  X  10      Hz 
0.  88  V 


=    9.  67  X  10     Hz/V 


f(Hz) 

f2 
(Hz2) 

B 
(Hz) 

V 

rms 
(MV) 

6v 

rms 

(Hz) 

(6i/rms)2 
(Hz2) 

S      (f) 
(Hz) 

60 
(dB)* 

5  X  103 

2.5  X  107 
(74  dB)t 

100 

660 

6.  38 

40.  7 

0.41 
(-3.9  dB)J 

-77.9 

2.5  X  103 

6.  2  X  106 

(68  dB) 

100 

7  80 

7.54 

56.9 

0.  57 
(-2.4  dB) 

-70.4 

1  X  103 

1.0  X  10 
(60  dB) 

100 

1000 

9.  67 

93.  5 

0.94 
(-0.  3  dB) 

-60.  3 

640 

4.  1  X  105 
(56  dB) 

10 

1200 

11.  6 

135.  0 

1.4 
(1.5  dB) 

-54.5 

320 

1.0  X  io5 
(50  dB) 

10 

4  60 

4.45 

19.  8 

2.  0 
(3.0  dB) 

-47.0 

210 

4.4  x  104 
(46  dB) 

10 

5  80 

5.  61 

31.5 

3.  2 

(5.  1   dB) 

-40.9 

150 

2.  2  X  104 

(44  dB) 

10 

680 

6.58 

43.  3 

4.  3 
(6.3  dB) 

-37.  7 

90 

8.  1  X  103 
(39  dB) 

1 

230 

2.  22 

4.93 

4.9 
(6.9  dB) 

-32.  1 

40 

1.  6  X  103 
(32  dB) 

1 

300 

2.90 

8.41 

8.4 
(9.  3  dB) 

-22.7 

20 

400 
(26  dB) 

1 

450 

4.  35 

18.9 

19.0 

(12. 8  dB) 

-13.2 

10 

100 
(20  dB) 

1 

700 

6.77 

45.  8 

46.  0 

(16.  6  dB) 

-    3.4 

*S.  .(f)   is  tabulated  in  deci 

°  2 

t  dB  relative  to   1  Hz 

J  dB  relative  to   1  Hz2  Hz"1 


2         -1 
bels  relative  to   1   radian     Hz 
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4.5.     Additional  Techniques  for  Frequency  Stability  Measurements 

at  X-Band 

It  has  been  found  convenient  and  desirable,   under  certain 
circumstances,    to  use  other  techniques  for  measuring  frequency  stabil- 
ity at  X-band.     Where  two  X-band  sources  are  available,   phase-  or 
frequency-locking  techniques  similar  to  those  used  at  HF  can  be  used. 
(See  figs.    15,    16,    and  17.  )    Good  wide-band  double -balanced  mixers 
with  coaxial  connectors  are  available  [11]  which  permit  many  of  the 
measurements  to  be  performed  -without  use  of  waveguide  components. 

The  measurement  setup  as  shown  in  the  block  diagram  of  figure  9 
can  also  be  used  at  microwave  frequencies  utilizing  a  computing  counter 
for  time  domain  measurements.     Extensive  measurements  of  frequency 
stability  have  been  made  on  stabilized  X-band  sources   [12].     Time 
domain  data  obtained  via  the  computing  counter  have  been  compared 
with  frequency  domain  data  obtained  via  several  methods. 

An  example  shown  in  Appendix  C  translates  the  frequency  domain 
data  of  figure  14  into  estimated  time  domain  performance  shown  in 
figure   18.     In  the  same  figure  we  have  plotted  time  domain  data  taken 
directly  via  a  computing  counter   (see  fig.    9). 

5 .     SUMMARY 

Terminology  and  concise  definitions  for  the  specification  of 
frequency  stability  have  been  given.     Recommended  techniques  and 
measurement  systems  were  described  in  detail  for  both  high  frequency 
and  X-band  signals.     Experimental  results  are  compared  using  various 
systems.     Examples  of  computations  are  also  included. 

A.    E.    Wainwright,    Howard  E.    Bell,    and  David  W.    Allan  have 
assisted  in  the  development  of  the  measurement  systems  reported  here. 
Their  contributions  are  gratefully  acknowledged.      The  authors  appreci- 
ate the  assistance  of  Mrs.    E.    Helfrich  in  the  preparation  of  the  manu- 
script. 
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APPENDIX  A 
Glossary  of  Symbols 

A  Peak-to-peak  voltage  of  a  beat  frequency  at  output  of  mixer 

ptp 

B  High  frequency  cutoff   f     (bandwidth) 

B  Analysis  bandwidth  (frequency  window)  of  the  spectrum 

analyzer 

f  Fourier  frequency  of  fluctuations 

f  Defined  as    B,    high  frequency  cutoff  (bandwidth) 

h  Coefficient  of   f        in  spectral  density  representation 

a. 

K  Calibration  factor  used  in  the  single  oscillator  stability 

measurement  system  for  microwave  frequencies, 

K     =     (Av)  /V 

rms       rms 

k,  n  Integers   (used  as  index  of  summation) 

s£(f)  Frequency  domain  measure  of  phase  fluctuations; 

Script    =£(f)    is  defined  as  the  ratio  of 

Power  density  (one  phase  modulation  sideband) 
Power   (total  signal) 

For  small    60,     S    „  (f)    m     2=?(f) 
60 

M  Total  number  of  data  values  available  (usually  M  »N) 

N  Number  of  data  values  used  in  obtaining  a  sample  variance 

P         ,  Total  power  of  signal 

total  F  e 

r  Parameter  related  to  dead  time;    r  =   T/T 

S       (f)  Spectral  density  of  frequency  fluctuations 

S       (f)  Spectral  density  of  voltage  fluctuations 

S  _  _  (f )  Spectral  density  of  phase  fluctuations;  S  _  .(f)   =  — =  S    (f) — — 

60  60  f2  y      f2 
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APPENDIX  A  cont. 

S  /  (l/1)  Spectral  density  of  the   (square  root  of  the)  radio  frequency 

V  RFP 

power 

5  (f)  Spectral  density  of   y  (Spectral  density  of  fractional  fre- 

quency fluctuations) 

T  Time  interval  between  the  beginnings  of  two  successive 

measurements 

t  Time  variable 

v  Root-mean- square  (noise)  voltage  at  output  of  mixer  as 

measured  by  a  spectrum  analyzer 

V  Root-mean- square  voltage  of  the  output  of  an  FM 

r  ms 

demodulator  due  to  intentional  modulation 

W  The  3-dB  resonance  linewidth  of  the  loaded  discriminator 

c 

cavity 

dx 
x  Time  interval  fluctuations;   — —    =    y,    hence    x  =   6T 

dt 

y  Fractional  frequency  fluctuations,     y    =     — 

o 
y  Average  of   y    over  a  specified  time  interval    T 

\    /  Time  average  operator  (usually  over  a  large  but  finite 

time  interval    T  ) 

A  Difference  operator 

6  Fluctuation  operator 
§V  Frequency  fluctuations 
60  Phase  fluctuations 

6V,  v  Voltage  fluctuations 

V  Signal  frequency  (carrier  frequency)  variable 

V  Average  frequency  of  source  (nominal  frequency) 
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APPENDIX  A  cont. 

cr  Square  root  of  a  variance 

2 
CT    (t)  Specific  Allan  variance  -where    N  =  2,    T  =   T 

y 

2 

cr    (N,  T,T,  f   )     Sample  variance  of  N  averages  of   y(t),  each  of  duration  T 

and  repeated  every   T    units  of  time   (Allan  variance) 

measured  in  a  post- detection  noise  bandwidth  of  £, 

h 

T  Sampling  time  interval 

T  Post-detection  averaging  time  of  the  spectrum  analyzer 

3, 

Cl  Signal  angular  frequency  (carrier  angular  frequency), 

CI    s    2-rrv 

CO  Fourier  angular  frequency  of  fluctuations,    CO   =    2tt£ 
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APPENDIX  C 

Translation  of  Data  from  Frequency  Domain 
Into  Time  Domain  Performance  Using  the  Conversion  Chart 

Referring  to  the  frequency  domain  plot  in  figure   14  it  is 

-3 
determined  that    S     ,  (f )    indicates  approximate    f        behavior  over  the 

50  _! 

total  range  plotted.     Therefore    S      (f)    is  nearly  proportional  to    f 

(i.e.  ,    flicker  frequency  noise).     At   f  =  1000  Hz,    S      (f)    is  equal  to 

-0.3  dB  relative  to   1  Hz   (see  Table   1).     The  carrier  frequency,    v   ,  is 

9.5  GHz. 

s    W    .     fffi£    =     "°"°-°3"*>     =     ,.04  X  lO"20^-1.  (CI) 

y  v  (9.5  X  10  Hz) 


o 


for   f    =    1000  Hz. 

(f) 

y 


S    (f)  — —      (see  conversion  chart)  (C2) 


h  =      S    (f)  Xf     =  (1.04  X  10_20Hz"1)  X  (103Hz)   =   10.  4  X  10-1' 

-1  y 

(C3) 
CT.2(T)       =      h       •    21n2     =     10.4  X  10"18  X  1.  39   =   14.5  X  10"18 


y 


(C4) 


O    (T)       =      3.  8  X  10"9   .  (C5) 

For  a  flicker  frequency  noise  there  is  no    T  dependence.     A  dashed  line 
at  this  calculated  value  is  plotted  on  the  same  graph  as  data  taken  directly 
in  the  time  domain  (fig.    18). 
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APPENDIX  D 
Spectral  Densities:     Frequency  Domain  Measures  of  Stability 

Donald  Halford 
Stabilities  in  the  frequency  domain  are  commonly  specified  as 
spectral  densities.     We  have  used  the  concept  of  spectral  density  ex- 
tensively in  the  preparation  of  this  report.      The  spectral  density 
concept  is  simple  and  very  useful,    but  care  must  be  exercised  in  its 
use.      There  are  at  least  four  different,    but  related,    types  of  spectral 
densities  which  are  used  in  this  paper.     In  this  Appendix,    -we  state  and 
explain  some  of  the  simple,    often-needed  relations  among  these  four 
often-used  types  of  spectral  densities. 

The  four  types  'which  we  have  used,    and  -which  are  most  relevant 
to  frequency  and  phase  fluctuations,    are 

S    (f)  Spectral  density  of  fractional  frequency  fluctuations 

(noise,    instability,    modulation).      The  dimensionality 
is  Hz      .      The  range  of   f   is  from  zero  to  infinity. 

S      (f)         Spectral  density  of  frequency  fluctuations   (noise, 

instability,    modulation).      The  dimensionality  is  Hz 
per  Hz.      The  range  of   f   is  from  zero  to  infinity. 

S     ,(f)         Spectral  density  of  phase  fluctuations   (noise,    insta- 
60  2 

bility,    modulation).      The  dimensionality  is  rad 
per  Hz.      The  range  of   f   is  from  zero  to  infinity. 

a?(f)  Script    =£(f)    is  a  frequency  domain  measure  of  phase 

fluctuations   (noise,    instability,    modulation).     Script 
s£(f)    is  defined  as  the  ratio  of  the  power  in  one  phase 
noise  sideband,    referred  to  the  input  carrier  fre- 
quency,   on  a  per  hertz  of  bandwidth  spectral  density 
basis,   to  the  total  signal  power,    at  Fourier  frequency 
f   from  the  signal's  average  frequency    v   ,    per  one 
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APPENDIX  D  cont. 

device.      The  dimensionality  is  Hz       .      The  range 
of  f   is  from  minus    v     to  plus  infinity. 

Each  of  these  spectral  densities  is  one-sided  and  is  on  a.  per 
hertz  of  bandwi dth  density  basis.      This  means  that  the  total  mean- 
square  fluctuation  (the  total  variance)  of  frequency,    for  example,    is 
given  mathematically  by 

co 

/S      (f )  df  , 
o 

and,    as  another  example,    since  Script    a?(f)    is  a  normalized  density, 
that 

+co 


/ 


^(f)  df 
V 

o 
is  equal  to  unity. 

Two-sided  spectral  densities  are  defined  such  that  the  range  of 
integration  is  from  minus  infinity  to  plus  infinity.      For  specification  of 
noise  as  treated  in  this  paper,    our  one-sided  spectral  density  is  twice 
as  large  as  the  corresponding  two-sided  spectral  density.      That  is, 

+co  +00  +00 

y*[sTwo-Sided]  df  =    zj  [sTwo-Sided]  ^       /^One- Sided  ]  ^ 

-co  o  o 

(Dl) 

Two-sided  spectral  densities  are  useful  mainly  in  pure  mathematical 
analysis.     We  recommend  and  use  one-sided  spectral  densities  for 
experimental  'work.     References   [1]  and  [2]  also  use  one-sided  spec- 
tral densities. 
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APPENDIX  D  cont. 


We  use  the  definition 


r-¥-.  m> 


V 
o 


and  it  follows  that 


V 

o 


(tj 


S    (f)     m      S£     (f)     =     /  — )      Sc     (f)  .  (D3) 


To  relate  frequency,  angular  frequency,  and  phase  -we  use 


2ir[v(t)3  =  0(t)  =  ^&   .  (D4) 


This  may  be  regarded  as  a  definition  of  instantaneous  frequency    v{t), 
From  equation  (D4),    a  direct  result  of  transform  theory  is 

„2  .    .2 


(i)  s6n(£)  ■  (?)  s6,<£)-  (D5: 


S6/»         '" 


Script  s£(f)  can  be  related  in  a  simple  way  to  S  (f),  but  only  for 
the  condition  that  the  phase  fluctuations  occurring  at  rates  f  and  faster 
are  small  compared  to  one  radian.      Otherwise  Bessel  function  algebra 

must  be  used  to  relate  Script  a?(f)    to    S      (f).      Fortunately,    the  "small 

6© 

angle  condition"  is  often  met  in  random  noise  problems.      Specifically 
we  use 


\2  rad2/ 


—  ]S60(f)   ,  (D6) 
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APPENDIX  D  cont. 

provided  that 

+  00 

/       S       ff')df    «    1  rad2  .  (D7) 

J  60 

For  the  types  of  signals  under  discussion  and  for     |f  |    <    V   ,    we 
use  as  a  good  approximation 

=£(-f)    »     =£(f)  .  (D8) 

Script    3?(f)    is  the  normalized  version  of   S  /  (V) ,    with  its  frequency 

parameter    f    referenced  to  the  signal's  average  frequency    v      as  the 
origin  such  that   f    equals    v   -    v   . 


Some  Mathematics  of  Phase  Sideband  Power  as  Related  to  Phase 

Fluctuations:     A  simple  derivation  of  equation  (D6)  is  possible.      We 

combine  the  derivation  with  an  example  which  illustrates  the  operation 

of  a  double-balanced  mixer  as  a  phase  detector.      Consider  two 

sinusoidal  5 -MHz  signals   (having  negligible  amplitude  modulation) 

feeding  the  two  input  ports  of  a  double-balanced  mixer.      When  the  two 

signals  are  slightly  out  of  zero  beat,    a  slow  sinusoidal  beat  with  a 

period  of  several  seconds  at  the  output  of  the  mixer  is  measured  to 

have  a  peak-to-peak  swing  of   A 

8  Ptp 
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APPENDIX  D  cont. 

Without  changing  their  amplitudes,    the  two  signals  are  retuned 
to  be  at  zero  beat  and  in  phase  quadrature   (that  is,     7T/2    out  of  phase 
with  each  other),    and  the  output  of  the  mixer  is  a  small  fluctuating  volt- 
age centered  on  zero  volts.      Provided  this  fluctuating  voltage  is   small 
compared  to    A   ,     /2,    the  phase  quadrature  condition  is  being  closely- 
maintained,    and  the  "small  angle  condition"  is  being  met. 

Phase  fluctuations 
,    respectively,    -where 


Phase  fluctuations    60  bet-ween  the  two  signals  of  phases    0      and 


T 


60    =    6(02  -   0j)  ,  (D9) 


will  give  rise  to  voltage  fluctuations    §V    at  the  output  of  the  mixer 

A 


6V   «      -$£  60,  (D10) 


where  we  have  used  radian  measure  for  phase  angles,  and  we  have 
used 

sin  60  m  60  (Dll) 

for  small    60   (60  <<  1  rad).     We  solve  equation  (D10)  for    60.    square 
both  sides  and  take  a  time  average 


<AptP) 


If  we  interpret  the  mean- square  fluctuations  of    60  and  of    6V", 
respectively,  in  equation  (D1Z)  in  a  spectral  density  fashion,    we  may 
■write 

44 


APPENDIX  D  cont. 


S6V(f) 

60  2 (A  )* 

rms 


where  we  have  used 


(A       )Z     =     8 (A  )2   ,  (D14) 

ptp  rms 


■which  is  valid  for  the  sinusoidal  beat  signal. 

For  the  types  of  signals  under  consideration,    by  definition  the 
two  phase  noise  sidebands   (lower  sideband  and  upper  sideband,    at    -f 
and   +f  from    v  ,    respectively)  of  a  signal  are  coherent  with  each  other. 
As  already  expressed  in  equation  (D8),    they  are  of  equal  intensity  also. 
The  operation  of  the  mixer  when  it  is  driven  at  quadrature  is  such  that 
the  amplitudes  of  the  two  phase  sidebands  add  linearly  in  the  output  of 
the  mixer,    resulting  in  four  times  as  much  power  in  the  output  as 
would  be  present  if  only  one  of  the  phase  sidebands  -were  allowed  to 
contribute  to  the  output  of  the  mixer.     Hence  for      f  I  <   V     we  obtain 


s6y<lfl>  5v^f1^o  +  £) 

— r  **       4  p  '  (D15) 

(A  )  total 

rms 


and,    using  the  definition  of  Script    sf(f), 


Vm^'o  +  £) 


*«      "      ^p— ~    ^S      (|£|)    ,  (D16) 

total 


provided  the  phase  quadrature  condition  is  approximately  valid. 
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APPENDIX  D  cont. 

The  phase  quadrature  condition  -will  be  met  for  a  time  interval  at  least 
T    long,    provided 


oo 

/ 


S       (f  ')df    <<     1  rad2,  (D17) 

(Znr)'1 


and  hence  equation  (D16)  is  useful  for  values  of   f   at  least  as  low  as 
(2ttT)       .      Equatior 
(D7)  respectively. 


(2ttT)       .      Equations   (D16)  and  (D17)  correspond  to  equations   (D6)  and 
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APPENDIX  E 

A  Sample  Calculation  of  Script  =£ 

Script    =£(f)    is  easily  measured  using  a  double-balanced  mixer 
as  a  phase  sensitive  detector,    together  with  a  spectrum  analyzer  'which 
can  measure  at  the  frequency   f.      Equation  (El)  is  valid  for  the  case 
•where  the  reference  signal  has  negligible  phase  noise  compared  to  the 
test  signal. 

_     2    /     one  unit  \     .  (El) 

n     \       ptp         / 

However,    eq  (E2)  is  the  valid  equation  when  we  have  two  equally  noisy 
signals   (test  and  reference)  driving  the  mixer. 


(E2) 


In  case  the  device  being  measured  has  amplification  of  times    n 

(i.e.  ,    frequency  multiplication  or  frequency  synthesis),    the  definition 

2 
of  Script    :£(f)    requires  that  the  factor    n      appear  in  eqs   (El)  and  (E2). 

For  example,    a  frequency  synthesizer  with  output  at  45.55  MHz,    and 

with  its  input  driven  at  5  MHz,    is  characterized  by   n    equal  9.11 

corresponding  to  a  phase  fluctuations  increase  expressed  in  decibels  of 

19. 2  dB. 


2 

10  log  n      =    20  log   n    =  20  log  9.11    =    19.  2  dB.  (E3) 

Both  (El)  and  (E2)  require  that  the  beat  signal  out  of  the  mixer  be 

sinusoidal  when  the  two  input  signals  are  slightly  out  of  zero  beat.      The 

peak-to-peak  amplitude  of  the  beat  is  A 

ptp 
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For  convenience  of  computation  and  plotting  it  often  is  advanta- 
geous to  set  the  beat  frequency  voltage  (before  locking)  to  some  special 

voltage  such  as  — , volts   (0.  316  V)  peak-to-peak  at  the  mixer  output. 

Then  (after  lock)  with  the  output  of  the  phase  detector  expressed  in  rms 
nanovolts  per  root  hertz,    direct  plotting  is  facilitated  for  Script    ;£(f)  in 
decibels  versus  frequency  in  hertz.     In  this  case   1000,    100,    and  10 
nanovolts  per  root  hertz  correspond  to  -110,    -130,    and  -150  dB 
respectively.     A  sample  calculation  demonstrating  this  convenience  is 
shown  below. 


Given: 


A,     =    0.316V  (i.e.  ,        1_    V),  (E4) 

ptp  ,/TcT 


v=     100  nV  •  Hz"1'2  rms  (E5) 


at   f    equal  to  20  Hz  for  a  pair  of  equally  noisy  devices  having  no  frequency 
multiplication  (n    =   1): 

*(f)  =  fVtwAounits  j 

V    APtp 

/  100  nVHz    /2\ 


0.  316  V 


yT^1 


(E7) 


lO'^Hz"1  =     -130  dB  (E8) 


at   f    equal  to  20  Hz.      Or,    using  logarithms; 
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„,r,        _„  ,            /     two  units  \ 
4{i)  =    20  log1()     — (E9) 

\         Ptp  / 


(10      V-  Hz    /2) 

=    20  log        * rjT =  20(-7  +  0.5)    =     -130  dB  (E10) 

10        (10    /a  V) 


at   f    equal  to    20  Hz. 

If  the  phase  noise  were  to  follow  flicker  law  (1/f),  at  f  equal  to 
1  Hz  the  mean  square  noise  would  be  20  times  worse  (13  dB  greater). 
That  is 

=£(1  Hz)    =  -130  dB  +  13  dB   =    -117  dB.  (Ell) 
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APPENDIX  F 


2 
A  Sample  Calculation  of  Allan  Variance,    (j    (T) 

— y — 


k=l 

(Fl) 


in  the  example  below: 

Number  of  data  values  available,   M  =  9 
Number  of  differences  averaged,    M  -  1   = 
Sampling  time  interval    T  =   Is 


Data  values 

First   differences 

(y) 

(^k+l   ~\] 

892 

809 

-    83 

823 

14 

798 

-    25 

671 

-127 

644 

-    27 

883 

239 

903 

20 

677 

-226 

M-l 

I 

k=l 


Fir 

st  differences 

squared 

6889 

196 

625 

16129 

729 

57121 

400 

51076 

(yk+1-y/     =        133165 


Based  on  these  data: 


a2{T)      =     ^T^    =     8322.81,  (F2) 

y  2(8) 


[0;    (T)]      =       x/8322.  81     =    91.  23,    N  =  2,    T  =  T  =  1  s.  (F3) 

In  this  example,    the  data  values  may  be  understood  to  be  expressed  in 

12 
parts  in  10 
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Using  the  same  data  as  in  the  above  example  it  is  possible  to 
calculate  the  Allan  variance  for    T   =  2  s    by  averaging  pairs  of  adjacent 
data  values  and  using  these  averaged  values  as  new  data  values  to 
proceed  with  the  calculation  as  before.     Allan  variance  values  may  be 
obtained  for    T  =    3  s    by  averaging  three  adjacent  data  values  in  a  similar 
manner,    etc.  ,    for  larger  values  of    T. 

Ideally  the  calculation  is  done  via  a  computer  and  a  large 
number,    M,    of  data  values  should  be  used.      (Typically  M  =  25  6  data 
values  are  used  in  the  NBS  computer  program.  )     The  statistical  confi- 
dence of  the  calculated  Allan  variance  improves  nominally  as  the   square 
root  of  the  number,    M,    of  data  values  used.      For    M  =  256,    the  confi- 
dence of  the  Allan  variance  is  not  expected  to  be  better  than  approxi- 
mately     ■    ,  X  100%   ~     7%      in  the   rms  sense.      The  use  of  M  »  1 
v/256 

is  logically  similar  to  the  use  of    B     •   T     »  1    in  spectrum  analysis 

a        a 

measurements,    where    B      is  the  analysis  bandwidth  (frequency  window) 
of  the  spectrum  analyzer,    and    T      is  the  post-detection  averaging  time 
of  the   spectrum  analyzer. 
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COMPUTING  COUNTER  a   (T)  PROGRAM 

y 

USING  AN  EFFICIENT  OVERLAPPING  ESTIMATOR 
David  W.    Allan 


(1) 
(2) 
(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 
(10) 
(ID 
(12) 
(13) 
(14) 
(15) 

(16)  bx'y 

(17)  +  (add) 


MANUAL 

(19) 

REPEAT 

Enter  carrier  or  basic  frequency 

(20) 

X  FER  PROGRAM 

(21) 

Nx  y 

[skip  to   (33)  if  program  is  already  in] 

(22 

N  x  y 

LEARN 

(23 

+   (add) 

CLEAR   x 

b^x 

* — 

(24 
(25) 

--* 

a      x 
«_• 

bx  y 

MODULE  or  PLUG-IN 
a      x 

(26 
(27 

axy 

-r   (divide) 

X  FER  PROGRAM 

(28 

^ 

MODULE  or  PLUG-IN 

(29 

c  x  y 

axy 

(30 
(31 

i      -=-   (divide) 
DISPLAY   x 

-   (subtract) 

(32 

PAUSE 

xy 

(33 

RUN 

X  (multiply) 

(34 

START 

(18) 


Vx 


Program  will  automatically  re- 
peat unless  righthand  PAUSE 
switch  is  in  HALT  position 


T       =        Sample  time   (computing  counter  "measurement  time") 
T  -   r     -    0.003  seconds  (compute  +  cycle  time) 

N     =    2 

Number  set  on  repeat  loop  corresponds  to  the  number  of  estimates  of  the 
variance.   For  good  confidence  levels   100  or  more  estimates  usually  are 
required. 
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Some  Important  References 

for 

Measurement  and  Specification  of  Frequency  Stability 

General  References 


1.  November  or  December  of  even-numbered  years  IEEE 
Transactions  on  Instrumentation  and  Measurement  (Conference  on 
Precision  Electromagnetic  Measurements,    held  every  two  years). 

2.  February  1966  Proceedings  of  the  IEEE,    Special  Issue  on 
Frequency  Stability  (IEEE-NASA  Symposium).    Very  useful. 

3.  Proceedings  of  the  IEEE-NASA  Symposium  on  the  Definition  and 
Measurement  of  Short-Term  Frequency  Stability  at  Goddard  Space 
Flight  Center,    Greenbelt,    Maryland,    November  23-24,    1964. 
Prepared  by  Goddard  Space  Flight  Center   (Scientific   &  Technical 
Information  Division,    National  Aeronautics  and  Space  Administra- 
tion,   Washington,    D.  C.  ,    1965).      Copies  available  for  $1.75  from 
the  U.    S.    Government  Printing  Office,    Washington,    D.  C.      20402. 

4.  The  annual  Proceedings  of  the  Symposium  on  Frequency  Control 
(Fort  Monmouth).      The  Proceedings  are  not  edited  nor  reviewed. 

5.  J.    A.    Barnes,    A.    R.    Chi,    L.    S.    Cutler,    et  al.  ,    "Characterization 
of  Frequency  Stability,  "  NBS  Technical  Note  394  (October   1970); 
also  published  in  IEEE  Trans,    on  Instr.    and  Meas.   IM-20,    No.    2, 
pp.    105-120   (May  1971).      This  is  the  most  definitive  discussion  to 
date  of  the  characterization  and  measurement  of  frequency  stability. 
It  was  prepared  by  the  Subcommittee  on  Frequency  Stability  of  the 
Institute  of  Electrical  and  Electronic  Engineers. 

6.  September   1968  IEEE  Transactions  on  Microwave  Theory  and 
Techniques,    Special  Issue  on  Noise. 

7.  Frequency  and  Time  Standards,    Application  Note  52,    Hewlett- 
Packard  Company  (November   19  65). 
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Some  Specific  Papers 

8.  D.    W.   Allan,    "Statistics  of  Atomic  Frequency  Standards,  " 
Proc.    IEEE,    vol.    54,    pp.    221-230,    February  1966. 

A  thorough  understanding  of  this  paper  is  important  for 
everyone  -who  •wishes  to  measure  and  quote  performance  of 
frequency  standards  in  the  time  domain,    e.g.,    a    versus    T 
plots.     The  data  analysis  must  take  into  account  the  number 
of  samples  taken  and  how  they  are  used. 

9.  David  W.   Allan,    B.    E.    Blair,    D.    D.    Davis,    and  H.    E. 
Machlan,    "Precision  and  Accuracy  of  Remote  Synchroniza- 
tion via  Portable  Clocks,    Loran  C,    and  Network  Television 
Broadcasts,  "  Proc.    25th  Annual  Symposium  on  Frequency 
Control,    Fort  Monmouth,    N.  J.  ,    pp.    195-208,    April  1971. 

10.  J.    R.    Ashley,    C.    B.    Searles,    and  F.    M.    Palka,    "The 
Measurement  of  Oscillator  Noise  at  Microwave  Frequencies,  " 
IEEE  Trans,    on  Microwave  Theory  and  Techniques  MTT-  1  6, 
No.    9,    pp.    753-760,    September   1968. 

11.  E.    J.    Baghdady,    R.    N.    Lincoln,    and  B.    D.    Nelin,    "Short- 
Term  Frequency  Stability:    Characterization,    Theory,    and 
Measurement,"  Proc.    IEEE,    vol.    53,    pp.    704-722,    July  1965. 
Among  many  other  topics,   the  possible  problem  of  AM  noise 
is  discussed. 

12.  Helmut  Brandenberger,    Frederic  Hadorn,    Donald  Halford, 
and  John  H.    Shoaf,    "High  Quality  Quartz  Crystal  Oscillators: 
Frequency  Domain  and  Time  Domain  Stability,  "  Proc.    25th 
Annual  Symposium  on  Frequency  Control,    Fort  Monmouth, 
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